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Due to the spin-orbit coupling (SOC) an electric current flowing in a normal metal or semi¬ 
conductor can induce a bulk magnetic moment. This effect is known as the Edelstein (EE) or 
magneto-electric effect. Similarly, in a bulk superconductor a phase gradient may create a finite 
spin density. The inverse effect, also known as the spin-galvanic effect, corresponds to the creation 
of a supercurrent by an equilibrium spin polarization. Here, by exploiting the analogy between a 
linear-in-momentum SOC and a background SU(2) gauge field, we develop a quasiclassical transport 
theory to deal with magneto-electric effects in superconducting structures. For bulk superconductors 
this approach allows us to easily reproduce and generalize a number of previously known results. 
For Josephson junctions we establish a direct connection between the inverse EE and the appear¬ 
ance of an anomalous phase-shift in the current-phase relation. In particular we show that ipo 
is proportional to the equilibrium spin-current in the weak link. We also argue that our results 
are valid generically, beyond the particular case of linear-in-momentum SOC. The magneto-electric 
effects discussed in this study may find applications in the emerging field of coherent spintronics 
with superconductors. 

PACS numbers: 74.50.+r Tunneling phenomena; Josephson effects - 74.78.Na Mesoscopic and nanoscale 
systems - 85.25.Cp Josephson devices - 72.25.-b Spin polarized transport 
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I. INTRODUCTION 

Over the past decades, superconductor-ferromagnetic 
(S-F) structures have been studied extensively [1, 2]. The 
spatial oscillatory behavior of the superconducting con¬ 
densate induced in the ferromagnet leads to interesting 
effects as oscillations of the density-of-state in F/S [3-5] 
and F/S/F [6] structures, oscillations of the Josephson 
current in S/F/S Josephson junction [7-10], and oscil¬ 
lations of the critical temperature [1]. Moreover, in the 
case of multi-domain ferromagnets or artificial multilayer 
structures with inhomogeneous magnetization, the sin¬ 
glet Cooper pairs from a superconductor can be trans¬ 
form into long-range triplet pairs that may explain the 
long-range Josephson coupling observed in S/F/S struc¬ 
tures [11-25]. Triplet correlations also leads to a depen¬ 
dence of the critical current on the magnetic configura¬ 
tion of diverse S/F structures [26-33]. Such phenom¬ 
ena suggest interesting perspectives of exploiting triplet 
correlations for the emerging field of coherent super- 
spintronics [34, 35]. Also promising applications might be 
found by using superconducting materials in combination 
with ferromagnetic insulators that may act as spin-filters 
[36-38]. In particular several thermal effects related to 
these material combination have been studied in recent 
works [39-44]. 


All the above mentioned phenomena in S/F structures 
originate from the interaction between the superconduct¬ 
ing correlations and the exchange field of the ferromag¬ 
net. However it has recently been shown that spin-orbit 
coupling (SOC) in S/F structures will also lead to, for 
example, a long-range triplet component [45, 46] and pe¬ 
culiarities in the density of states [47-49]. On the other 
hand, transport properties of non-superconducting struc¬ 
tures with strong SOC are being intensively studied be¬ 
cause of their potential application in a novel direction 
of spintronics, which exploits the coupling between spin 
and charge currents [50-53]. 

In particular, the SOC in semiconductors and normal 
metals is at the root of a number of interesting phenom¬ 
ena that originate from the coupling between the charge 
and spin degrees of freedom. Prototype of these phenom¬ 
ena is the spin Hall effect (SHE) [54-64] which consists 
in the creation of a spin polarized current by an electric 
field. Reciprocally, by means of the inverse SHE a spin 
current can create an electric field [65-67]. These effects 
allow to generate and detect spin polarized currents in 
non-magnetic materials [68-72]. 

There is also another relevant effect in normal systems 
related to the SOC. It consists in creating a stationary 
spin density S a , along the a-direction in response to an 
electric field Ek applied in ^-direction [73, 74]. Within 
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linear response, this effect is described by 

S a (co) = a a k (co)E k (uj) , (1.1) 

where the sum over repeated indexes is implied here, and 
throughout this paper. In particular, in 2D systems with 
Rashba SOC, the applied electric field and the gener¬ 
ated spin density are perpendicular to each other. This 
magneto-electric effect, also called the Edelstein effect 
(EE), has been observed in experiments [60, 75]. The 
Edelstein conductivity a k (cj) in Eq.(l.l) is related to the 

Kubo correlator x k (aO = spin an d 

current operators via a k (cj) = xt (, u ) /' 1UJ [76]. Because 
of the gauge invariance in normal systems the function 
X k (cj) should vanish in the limit uj 0 reflecting the 
fact that there is no response to a static vector potential. 
Therefore the a k (0) = cr k remains finite and describes 
the dc EE. It has been pointed out in Ref. [76] that this 
property, together with the Onsager reciprocity princi¬ 
ple, implies that the inverse dc EE, also referred to as 
the spin-galvanic effect, consists in generating a charge 
current j k by a steady spin generation induced by a time- 
dependent magnetic field via the paramagnetic effect: 


Jk = cr k 


9HBB a 


( 1 . 2 ) 


with the Lande < 7 -factor, pb the Bohr magneton, and 
B a the time derivative of the magnetic field along the 
a-axis. The inverse EE effect has also been observed in 
experiments [77, 78]. 

Similar magneto-electric and spin-galvanic effects 
should also exist in superconductors [79, 80]. However, 
there the physical situation is different because in the 
presence of the superconducting condensate the gauge 
invariance does not forbid the existence of a finite static 
current-spin response function x%- I n con t ras t 1 ° the 
normal case, in a superconductor an equilibrium electric 
(super-)current can flow in the absence of an external 
electric field. The supercurrent j = n s v s (here n s is 
the density of superconducting electrons and v s the su¬ 
perfluid velocity) is proportional to the gradient of the 
macroscopic gauge-invariant phase V(p = \7(p — eA ^ v s , 
which is the physical field coupled to the current operator 
in the Hamiltonian of a superconductor. The existence 
of such a gauge-invariant field implies that the static 

response function xt = can nonzero 

without violating the gauge invariance. In principle, a 
supercurrent can thus generate an equilibrium spin po¬ 
larization according to the general linear response rela¬ 
tion: 


= X%dk<P , (1.3) 

where d k = d/dx k . This effect has been indeed theoret¬ 
ically demonstrated by Edelstein for a 2D superconduc¬ 
tor with Rashba SOC, who calculated the proportionality 
tensor xt a ^ temperatures T close to the critical super¬ 
conducting temperature T c , in both pure ballistic [79] 
and diffusive [80] limits. 


Because in the superconducting state the response 
function x% at cj 0 is finite, the reciprocity of the 
EE effect becomes complete. In contrast to the normal 
case, in superconductors a static Zeeman field B can in¬ 
duce a supercurrent j k . Therefore, instead of Eq. (1.2) 
the following relation holds 

jk = exih a , (1.4) 

where h a = An explicit expression of this type 

has been obtained in a particular case of a 2D ballistic 
superconductor with intrinsic Rashba SOC [81, 82]. 

It is then clear that the free energy of a superconductor 
with a SOC must have a term of the Lifshitz-type 

Fl = h a X a k d k (p , (1.5) 

and equations (1.3)-(1.4) follow directly from the gen¬ 
eral thermodynamic definitions of the spin and current 
densities, S = 5F/Sh and j = —SF/SA. 

In principle, equations (1.3) and (1.4) apply for bulk 
superconductors, but one can expect similar effects to 
occur also in a S-X-S Josephson junction, between two 
massive superconductors (S) and a normal or ferromag¬ 
netic bridge X with an intrinsic SOC. In a Josephson 
junction the supercurrent depends on the phase differ¬ 
ence cp between the superconducting electrodes. In the 
particular cases of a weak proximity effect between the S 
and the X, or in the high-temperature regime (T < T c ), 
the current phase relation is given by j = j c sin ip, where 
j c is the critical Josephson current. 

When the SOC competes with a Zeeman effect, the 
natural conjectures following Eqs. (1.3)-(1.4) are: (i) In 
accordance with Eq. (1.3), the flow of a supercurrent 
may generate a spin polarization in the X bridge (the 
Edelstein effect); (ii) In turn, from Eq. (1.4), a Zeeman 
(spin-splitting) field may induce a supercurrent through 
the junction, even if the phase difference between the 
electrodes vanishes (the inverse Edelstein effect). 

In other words, the inverse EE is presumably the cause 
of an anomalous phase ipo which modifies the current 
phase relation according to j = j c sin (cp — </?o)> with a 
non-trivial (i.e., non equal to 0 or n) equilibrium phase 
ipo- This defines the so called </?o-junctions, a subject 
that has been extensively studied in the past years in 
different systems, including conventional superconduc¬ 
tors with SOC [83-95], with triplet correlations [96-101] 
or in contact with topological materials [102, 103], and 
also in hybrid systems with non-conventional supercon¬ 
ductors [104-112], quantum dots [113-115], and hybrid 
(0 — 7 r)-structures [116-119]. (^o-junctions may produce 
a self-sustained flux when embedded in a SQUID ge¬ 
ometry [ 120 ], act as phase batteries in coherent circuits 
[ 121 , 122 ], present a current asymmetry and act as a su¬ 
percurrent rectifiers [113]. 

In the present work we develop a complete theory of 
the magneto-electric and spin-galvanic effects in hybrid 
superconducting structures and confirm the above con¬ 
jectures. We focus on systems with linear in momentum 




3 


SOC that can be conveniently described in terms of an 
effective background SU(2) gauge field. This allows us to 
use the SU(2) covariant quasiclassical equations for the 
Green’s functions derived in Ref. [45, 46, 123]. We estab¬ 
lish a connection between the tensor xt i n Eqs. (1.3-1.4) 
and the equilibrium spin-current 3°j [124, 125]. We show 
that in a generic S-X-S Josephson junction the condi¬ 
tion for a nontrivial anomalous phase <po to appear is 
that 7 ^ 0, where h a can be either an external Zee- 
man field or the internal exchange field of a ferromagnet. 
Our SU(2) covariant formulation results in a simple and 
tractable system of equations to describe hybrid struc¬ 
tures with arbitrary linear in momentum SOC, temper¬ 
atures, degree of disorder, and quality of the hybrid in¬ 
terfaces. We also show that qualitatively our results are 
generically valid beyond the particular case of the linear 
in momentum SOC. 

The structure of the paper is the following: In the next 
section we present a qualitative discussion of the super¬ 
conducting proximity effect in structures with SOC and 
its connection with the spin diffusion in normal systems. 
This qualitative analysis allows us to guess the form of 
the quasiclassical equations for superconducting struc¬ 
tures in the presence of generic spin fields, and in particu¬ 
lar to explicitly show the analogy between the charge-spin 
coupling in normal systems and the singlet-triplet cou¬ 
pling in superconducting ones. In section III we present 
our model, discuss the associated symmetries, and derive 
microscopically the quasiclassical equations for generic 
linear in momentum SOC. In section IV we use the de¬ 
rived equations to explore the magneto-electric effects 
in bulk superconductors. We generalize the previously 
known results for the EE and its inverse obtained for 2D 
Rashba SOC [79-81, 126] to generic linear-in-momentum 
SOC, and relate them to the spin current and the SU(2) 
gauge-fields. In section V we explore the Josephson ef¬ 
fect through a S-X-S diffusive junction and in section VI 
through a ballistic one. In both cases we show that the 
anomalous phase (po is proportional to 3^h a and deter¬ 
mine its dependence on other parameters of the struc¬ 
ture, like temperature and length. We finally present our 
conclusions and discuss possible experimental setups to 
verify our predictions in Section VII. 


II. DIFFUSION OF SUPERCONDUCTING 
CONDENSATE IN THE PRESENCE OF 
SPIN-ORBIT COUPLING: HEURISTIC 
ARGUMENTS 

Before presenting the full quantum kinetic theory it is 
instructive to discuss at the qualitative level the main fea¬ 
tures of the proximity induced superconductivity in the 
presence of an intrinsic SOC. For this sake we present 
a simple heuristic derivation of the equations describing 
the coupled motion of the singlet and triplet components 
induced in a ferromagnet from a bulk 5-wave supercon¬ 
ductor. 


Let us consider a S-X-S junction, where X is a diffusive 
ferromagnet. We assume that the system is at equilib¬ 
rium, and that the proximity effect between S and X is 
weak. In such a case the junction is fully described by the 
quasiclassical anomalous Green function /(r), which de¬ 
scribes the superconducting condensate in X. In general 
/(r) is a 2 x 2 matrix in the spin space / = f s l + / t a cr a . 
Here the scalar f s and the vector with components / t a de¬ 
scribe the singlet and the triplet components of the con¬ 
densate, respectively. In this section we show, that the 
functions f s ( r) and f t (r) are reminiscent of the charge 
and spin density in the normal systems. 

In the absence of SOC, but in the presence of the ex¬ 
change field h the diffusion of the condensate is described 
by the well known linearized Usadel equations (see e.g. 
Ref. [ 2 ]), 

DX/ 2 f s - 2 |w„| f s + 2i sgn(co n )h a ft = 0 , (2.1) 

-DV 2 /t a “ 2 \oj n | ff + 2 isgn(w n )/i°/ s = 0 , ( 2 . 2 ) 

where D is the diffusion constant, and uj n is the Mat- 
subara frequency. The terms proportional to 2|cj n | are 
responsible for the decay of the superconducting corre¬ 
lations in the normal metal. The last terms in the left 
hand sides of Eqs. (2.1)-(2.2) describe the usual singlet- 
triplet coupling coming from the exchange field. It is 
worth emphasizing the presence of imaginary unit i in 
the exchange field terms, which reflects the breaking of 
the time reversal symmetry. Because of this, the singlet- 
triplet conversion due to the exchange field is always ac¬ 
companied with a phase shift of 7 r/ 2 . This point will be 
of primary importance in the following for understanding 
the origin of the anomalous phase (po . 

To understand how the Usadel equations (2.1)-(2.2) are 
modified in the presence of SOC we recall the description 
of the diffusion of spin S(r) and charge n( r) densities in 
normal systems. The general spin diffusion equation in a 
normal conductor with SOC takes the form, 

d t S a - DV 2 S a = 7 a , (2.3) 

where T a is a so called spin torque. In the absence 
of SOC, 7 a = 0 and hence spin is a conserved quan¬ 
tity which satisfies the usual spin diffusion equation. In 
non-centrosymmetric materials SOC acts as an effective 
momentum-dependent Zeeman field that causes preces¬ 
sion of spins of moving electrons. This precession breaks 
conservation of the average spin, and shows up formally 
as a finite torque 7 a 7 ^ 0 in Eq. (2.3). In the diffusive 
regime the motion of the electrons consists of a random 
motion superimposed on an average drift caused by the 
density gradients. The spin precession related to these 
types of motion generate the corresponding contributions 
to the spin torque. To the lowest order in gradients the 
general expression for the torque can be written as fol¬ 
lows [127-129], 

7“ = D [-r ab 5 b + 2 Ppd k S b + C%d k n] . (2.4) 
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Here the first term describes the Dyakonov-Perel (DP) 
spin relaxation that originates from the spin precession 
of randomly moving electrons [54]. The positive definite 
matrix T ab is the DP relaxation tensor with the eigen¬ 
values equal to the inverse squares of the DP spin relax¬ 
ation lengths. The other two contributions to the torque 
are related to the average motion of spins. In particu¬ 
lar, the second term in the right hand side of Eq. (2.4) 
originates from the diffusive motion of spins caused by in¬ 
homogeneities of the spin density distribution. The cor¬ 
responding spin precession is described by antisymmetric 
(spin rotation) matrices Pg b = —P% a with ||P|| ~ l/£ so , 
where £ so is the spin precession length. 

The last term in Eq. (2.4), which is proportional to 
the charge density gradient, can be called the spin-Hall 
torque. The charge density gradient generates the charge 
current which is then transformed to the spin current via 
the spin Hall effect. Precession of the spins driven by 
the charge density gradient, via the spin Hall effect, is 
the origin of the spin-Hall torque in Eq. (2.4). The spin- 
Hall torque is parameterized by the tensor C% which is 
proportional to 0 s n/£ scn where # s h is the spin Hall angle 
- the conversion coefficient between the charge and the 
spin currents. 

Equation (2.3) with the spin torque of Eq. (2.4) is com¬ 
monly used in spintronics context to describe spin dy¬ 
namics in semiconductors with intrinsic SOC [127-129] 
(for a discussion between intrinsic and extrinsic SOC, see 
e.g. [64]). In the stationary case the diffusion equations 
for the spin and charge densities reduce to 


V 2 n + Cld k S a = 0, (2.5) 

V 2 S a - r 06 ^ 6 + 2 Pfd k S b + Cld k n = 0. (2.6) 

It is important to emphasize here that spin-charge cou¬ 
pling mediated by the spin-Hall torque (C%) is responsi¬ 
ble for the EE. This can be seen directly from Eq. (2.6): 
A uniform charge density gradient produces a uniform 
spin density given by S a = (T~ 1 ) ab C^.dkn. 

We can now construct the Usadel equations in the pres¬ 
ence of SOC in analogy to the normal case. Since SOC 
does not violate the time reversal symmetry it acts in ex¬ 
actly the same way on the time-reversal conjugated states 
composing the Cooper pair. Therefore the diffusion of the 
singlet and the triplet condensates should be modified by 
SOC in complete analogy with the diffusion of the charge 
and spin densities in normal systems. The formal con¬ 
nection between the diffusion of the triplet condensate 
function / t a in superconductors and the spin density S a 
in normal metals has been discussed recently in Ref. [46] , 
and it has been also noticed in Ref. [86]. Hence, in order 
to include the effects of SOC in the Usadel equations all 
we need to do is to replace the diffusion operators (the 
Laplacians) in Eqs. (2.1) and (2.2) with the diffusion op¬ 
erators entering Eqs. (2.5) and (2.6), respectively. The 
result is the following system of equations describing a 
coupled diffusion of the singlet and triplet condensates 


in the presence of SOC, 


V 2 /, 

V 2 /“ 


r 2 h a 

n 2 u f s + sgn (w n ) [i —/“ + C%d k f? 
(Kl6 ab +T ab )f! + 2P£ b d k f! 


r 2 h a 

sgn (u„) |^i—/ s + C k d k f s 


= 0 . 


= ((2.7) 

( 2 . 8 ) 


In contrast to the normal case, in addition to the DP 
relaxation, both the f s and ft experience an additional 
d ecay prop ortional to the inverse decay length = 
yj2\u) n \/D, due to the finite lifetime of the superconduct¬ 
ing condensate in the normal metal. 

The most important novel feature of Eqs. (2.7)-(2.8) 
is the presence of two mechanisms for the singlet-triplet 
coupling which are described by the two terms in the 
square brackets. The first mechanism is the above dis¬ 
cussed Zeeman coupling related to the modification of 
the internal structure of the Cooper pair by the spin¬ 
splitting field h [see Eqs. (2.1)-(2.2)]. The second chan¬ 
nel of singlet-triplet coupling comes from the spin-Hall 
torque, which converts the gradient of f s into f t and 
vise versa, in a complete analogy with the EE in normal 
systems. The corresponding singlet-triplet “conversion 
amplitudes” have a relative phase shift of 7r/2, which is 
related to the different transformation properties of the 
Zeeman and spin-orbit fields with respect to the time 
reversal. We will see in the next sections that the inter¬ 
ference of these two singlet-triplet conversion channels is 
indeed responsible for the magneto-electric/spin-galvanic 
effects in superconductors, and, in particular, for the ap¬ 
pearance of the intrinsic anomalous phase ipo in Joseph- 
son junctions. 

Although the present heuristic derivation of Eqs. (2.7)- 
(2.8) may seem imprecise, it uncovers a simple, but deep 
connection between the physics of inhomogeneous super¬ 
conductors with SOC and the well known spintronics ef¬ 
fects, such as the spin Hall effects and direct and inverse 
magneto-electric effects (EE). In Sec. Ill we present a 
rigorous derivation of the quasiclassical kinetic equations 
for superconductors with a linear in momentum SOC, 
which in the diffusive limit confirms the correctness of 
Eqs. (2.7)-(2.8). In the rest of the article we study in de¬ 
tail the physical consequences of the interference of the 
two singlet-triplet conversion channels and their connec¬ 
tion with the theory of (^o-Josephson junctions. 


III. THE MODEL AND BASIC EQUATIONS 


In this section we introduce our model and discuss 
the symmetries associated with superconducting systems 
in the presence of spin-orbit coupling (SOC). We also 
present the derivation of the quasiclassical equations in 
the presence of linear in momentum SOC. 
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A. The Hamiltonian in the presence of generic 
SOC and symmetry arguments for the appearance 
of an anomalous phase 


Our starting point is a general Hamiltonian describing 
a metal or a semiconductor with a linear in momentum 
SOC, an exchange field and superconducting correlations 


Hm dr 


j dr T ] , (3.1) 


where ( r ) are the annihilation operators for spin up 
and down at position r, and ^ is the spinor 

of creation operators. Hq is the free electron part [149] 

Ho= (P±^_Ai)__ Ao + V ( 3 . 2 ) 

2m 

where p the chemical potential and VJ mp the potential in¬ 
duced by non-magnetic impurities. The magnetic inter¬ 
actions appear in two places: as a SU(2) scalar potential 
Ao = Agcr a /2, describing for example the intrinsic ex¬ 
change field in a ferromagnet or a Zeeman field in a nor¬ 
mal metal, and as a SU(2) vector potential Ai = A^a a /2^ 
describing the SOC. The latter is associated to the mo¬ 
mentum operator[150] Pi = — idi in the form of a minimal 
coupling pi — Ai. In practice, all the linear-in-momentum 
SOC can be represented as a gauge potential (see e.g. 
[130] or [131] and references therein). In the widely 
studied case of a free electron gas with Rashba SOC, 
A v x = —Ay = — a , while in the presence of Dresselhaus 
SOC A x x = — A v y = /3. Finally, V = V (r) < 0 in the sec¬ 
ond term of the r.h.s of Eq. (3.1) describes the coupling 
strength which gives rise to superconductivity in some 
regions of space. 

In analogy to electrodynamics one can define the four- 
potential Ap, with space components (/i = 1, 2, 3 or p = 
x, y , z) given by the SOC and the time component (p = 0) 
by the Zeeman field. Following the analogy one can define 
the strength tensor 

3^ = l^o- 0 = d^Av - d- i[An,A v \, (3.3) 


and the electric and magnetic SU(2) fields 


= %k and B? = e ijk ?“ k , (3.4) 

where Eijk is the Levi-Civitta symbol. 

In normal metals and semiconductors, the SHE and 
EE are consequences of the existence of a finite SU(2) 
magnetic field. For a pure-gauge vector potential the 
SOC can be gauged out [46], the SU(2) magnetic field is 
zero, and hence the SHE and EE do not appear [151]. 
Following our analogy, in the superconducting case an 
anomalous phase can only appear if the SU(2) magnetic 
field is finite. This explains why S-F-S junction without 
SOC do not present any magneto-electric effect, or equiv¬ 
alently, no anomalous phase. As it is well known, the 


ground state of S-F-S junctions corresponds to a phase 
difference either equal to 0 or to 7r [1, 2, 132]. 

A simple way to describe qualitatively magneto-electric 
effects in a superconductor is to provide simple symmetry 
arguments. Let us consider a ballistic superconductor at 
T close to its critical temperature T c and focus on the 
Ginzburg-Landau free energy. In such an expansion, a 
SOC is responsible for the presence of a first-order deriva¬ 
tive of the order-parameter, the so-called Lifshitz invari¬ 
ant which describes most of the original phenomenology 
of non-centrosymmetric superconductors [82]. Assuming 
that the amplitude of the order parameter is constant but 
its phase position-dependent, the Lifshitz invariant reads 
Fl oc Tidip where Ti is a vector which has to be odd 
with respect to the time-reversal operation, and SU(2) 
invariant. As discussed in Ref. [92], to the lowest order 
in SOC the Lifshitz invariant for a superconductor can 
be expressed in terms of the SU(2) fields: 


F l oc Tr {T 0i %} dpp = (£ a x £ a ) • Vp . (3.5) 

If we focus on the static case, the electric field is given 
by 3^oj = — djAo. Moreover we define the equilibrium 
spin current [125] in terms of the SU(2) magnetic field as 
3j = = d^ij/dxi — i [Ai, $ij\. If Aq is spatially ho¬ 

mogenous, for example induced by an external magnetic 
field, Eq. (3.5) reads [92] 

F l oc Affidup . (3.6) 

This Lifshitz invariant agrees with the ones derived from 
microscopic considerations [133] or quasi-classic expan¬ 
sions [134] for a particular sort of SOC. 

Eq. (3.6) confirms our guessed Eq. (1.5) and demon¬ 
strates that the Edelstein response tensor xt behaves like 
the spin current tensor df- The form of Fl in Eq. (3.6), 
in terms of the equilibrium spin current, suggests that 
our results remain valid for any momentum dependence 
of the SOC. We now proceed to derive the quasiclassical 
equations and provide a microscopic description of the 
magneto-electric effects in superconductors. 


B. The quasiclassical equations in the presence of 
SOC 

In order to describe the transport properties of hy¬ 
brid structures containing superconducting, normal (N) 
and/or ferromagnetic (F) layers with interfaces, arbitrary 
temperature and degree of disorder we have to go beyond 
the Ginzburg-Landau limit. We present here the quasi¬ 
classical equations [135-138] for the Green’s functions in 
the presence of a non-Abelian gauge-field [45, 46, 123] 
(for a similar discussion in normal metal, see [139]). 

We follow here the derivation presented in Ref. [46]. 
The basic transport equation derived from Hamiltonian 
(3.1) for the Wigner-transformed covariant Green func- 
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tions G (p, r) in the time-independent limit reads : 


73 (<^n — h^o) — iA + 7^-, G 

ZT 

0’ ( 3 - 7 ) 

where u; n = 2T7r(n + l/2) is the fermionic Matsubara 

( 0 e i(p \ 

frequency, A = A ( ^ 1 is the (5-wave) gap pa¬ 

rameter of amplitude A and phase ip. The scattering at 
impurities is described within the Born approximation, 
where r is the elastic scattering time, (g) is the GF ma¬ 
trix integrated over the quasiparticle energy, and (• • •) 
describes the average over the Fermi momentum direc¬ 
tion. 

After integration of (3.7) over the quasiparticle energy 
and by using the fact that G is peaked at the Fermi level 
one obtains the generalized Eilenberger equation [46, 92]: 

v F («iVj) g + [73 (uj n - LA 0 ) - iA, g] 

= • (38) 

where n*, i = x,y,z are the components of the Fermi 
velocity vector. When deriving (3.8) we have neglected 
corrections to the exchange term Ao of the order of 
\Aj\/pp < 1. In fact, one sees from (3.7) that 
{73T0 j,dG/dpj} scales like {Ajd/dpj, — i [73/I0, G\ } 
since To j — — i [Ao,Aj], and so it renormalizes the term 
—i [r3Ao,G] already present in (3.7). The correction to 
Aq is of the order Aj/pp <C 1 and we neglect them from 
now on. 

In the Nambu space g reads 

*=(-/-»)' <39) 


-V,G + 

rri 


C. Linearized quasiclassical equations in diffusive 
and pure-ballistic limits 


In the present work we mainly consider two limiting 
cases: the pure ballistic one in which r —)> 00 and the 
diffusive limit where r is a small parameter. The trans¬ 
port equation in the ballistic limit is directly obtained 
from (3.8) by neglecting the right-hand side. The diffu¬ 
sive limit is a bit more puzzling. Because of the anti¬ 
commutator in the l.h.s of Eq. (3.8), the normalization 
condition g 2 = 1 does not hold directly and therefore 
the usual derivation of the Usadel equations cannot be 
carried out [140]. There is, however, a way out of this 
puzzle if one assumes that the amplitude of the anoma¬ 
lous GF’s, / in (3.9) is small. Then the matrix GF (3.9) 

0 / 

-/ 0 

linearized Eilenberger equation becomes an equation for 


can be written as g ~ sgn(cj n )r3 + 


and the 


( vpniVi + 2w ra ) / - {LAo, /} + 2iA sgn (w„) + 

= (/ - (/>) (312) 

This linearization procedure is justified in two cases: ei¬ 
ther for temperatures close to the critical temperature T c 
when the amplitude of the order parameter A is small, or 
in S-X structures when the proximity effect is weak due 
to a finite interface resistance for arbitrary temperature. 

In the diffusive limit one can expand / « / 0 + rikfk + 
• • •, in angular harmonics where (/) = fo fk- We first 
average (3.12) over the momentum direction: 

V J7 1 ~ 

-jr^Vfc/fe + - LA 0 , fo} = —2iA sgn (u>„) , (3.13) 

where dim = 1,2,3 is the dimension of the system. Next 
we multiply Eq.(3.12) by rik and average over the mo¬ 
mentum direction to obtain 


where the g, / are matrices in the spin space which de¬ 
pend on the spaces coordinates the momentum direc¬ 
tion rii and the Matsubara frequency. The time-reversal 
conjugate g and / are defined as g( n) = cr y g*(—n)a y and 
/ = cr y f*(—n)a y . The latter is the anomalous GF which 
describes the superconducting correlations. 

From the knowledge of g one can calculate the charge 
current (density) 


j 


ineNoT 

2 


X ^ ( vfT3 9 ) 

UJ n 


(3.10) 


with e the electron charge and Nq the normal density 
of states for each spin. Whereas the spin polarization is 
given by 


S = 


i7rA 0 T 


X Tr (t-3 erg) 


UJn 


(3.11) 


VF V/c/o 4~ iAq •) fk\ 2— 

(3.14, 

T 

Eqs. (3.13) and (3.14) constitute a closed set of coupled 
differential equations for /o and fk . In particular from 
Eq. (3.14) we can write fk in terms of fo up to terms of 
second order in r: 

f k « -rsgn(w n )ti F V fc /o - 7-2 77 \^kj, X/o} 

+ t 2 v f [w„ - LAo, Vfe/oj H-• (3.15) 

Note that the Usadel equation was obtained in several 
works in the absence of gauge-fields, where one skipped 
the terms of the order r 2 . We keep here these terms since 
they are crucial for the description of magneto-electric 
effects [86, 92, 141]. 


2 








7 


The equations can be further simplified by noticing 
that the anti-commutator in the second line of Eq. (3.15) 
can be written as 


jw n - LA 0 , V fe /o| = Vfe K - LA 0 , /o} 

+ i{v fc A 0 ,/ 0 } • (3.16) 

In virtue of (3.13), the first term in the right-hand-side 
of the last equation is in fact of order r and so this term 
in (3.15) is of order r 3 and can be neglected. The sec¬ 
ond term reads V/c^Ao = —i [-Afe, -Ao] = ^/co for a space- 
independent gauge-potential. This electric field renor¬ 
malizes the paramagnetic effects .Ao, and is neglected in 
the following. Finally, we replace (3.15) into (3.13) to 
obtain the Usadel equation for /q: 


- sgn (ui n ) £>V 2 /o + {uj n - LA 0 , /o} 

- 7 ^ {ViJjj, Vj/o} = — 2 iAsgn (w n ) (3.17) 

with D = dim the diffusion constant. This equation 
is supplemented by the generalized Kupriyanov-Lukichev 
boundary condition [142] 




Vi/o- 


rsgn (cj n ) 
2m 


{^V/o} 


= —7 ,/bcs 


(3.18) 

at an interface located at position xo between a bulk su¬ 
perconductor described by the anomalous GF /bcs and 
the X bridge. The interface is characterized by the trans¬ 
parency 7 and normal vector of component N*. For a 
fully transparent interface, we impose the continuity of 
the GFs. 

We now need to write the current and spin density in 
terms of the isotropic anomalous GFs. It is easy to verify, 
by checking its conservation, that in the linearized case 
the electric current, Eq. (3.10), is given by: 


3 = 17re ^° r ^ Tr (v F ff) s gn u n , (3.19) 

^n 


and correspondingly in the diffusive limit 


ji = ineNoDT ^ IV {/oV/o - AiV/o + 

UJ n 

+ TSS 2 n ^ ra) (/o {*a,Vifo} + f 0 Vi/o}) } • 

(3.20) 

The spin polarization (3.11) is more subtle to deal with 
in the linearized approximation, since the normalization 
condition do not apply in our case. In accordance with 
the case without SOC, one may assume that it can be 
expressed in terms of the isotropic anomalous / as: 

S a = inN 0 T Tr 7 °//} sgn (w„) (3.21) 


with (cy a ff) = cr a fofo in the diffusive limit. In the next 
section we will show a posteriori that these expressions 
leads to the known results in bulk systems in the presence 
of Rashba SOC. 

For the following discussions it is convenient to write 
the anomalous GF / as the sum of singlet (scalar) and 
triplet (vector in spin space) / = f s + / t a cr a , and to ex¬ 
pand all the spin variables in term of Pauli matrices: 
3^ij = 3^ijCr a / 2, = A%a*/2. From Eqs. (3.12) we ob¬ 

tain the equations for the singlet and triplet components 
in the ballistic case: 


( v F riidi + 2 u n ) f s 


— 2 isgn (w„) A 

+ 7 +^^)/.“ p' 22 ) 


f ~ \ a ( rqTf- d \ 

<>F», (V,/,) + 2w„/ ( = (MS + 

(3.23) 

Equivalently, from Eq. (3.17) one obtains the equa¬ 
tions for the isotropic part of the singlet f s o and triplet 
/to components in the diffusive case (for simplicity we 
skip the subindex 0 ) : 


(a? - 4) /. - 2 ig+ 


+ sgn (uj n ) 

(tv/" - 7/* a 


A a i- 

(3idi 

D 2m 


(\ a 


ft= 0 (3.24) 


+ sgn (w„) (i-^- + T^dfdi) f s = 0 (3.25) 


We write the covariant derivative as V* = di — i , ] = 
di + Pi , where Pi is a tensor dual to Ai with components 
pab _ £ abcj^C' Thus, ViVi = d\ + 2 — f, where 

PiPi = —f. By noticing that (r/2m)3^ = Cjf, it is easy 
to verify that the diffusive equations Eqs. (3.24-3.25) 
are identical to those derived in section II from heuristic 
arguments [Eqs. (2. 7-2. 8 )]. One should emphasize though 
that while Eqs. (3.24-3.25) are derived for the particular 
case of linear in momentum SOC, Eqs. ( 2 . 7-2. 8 ) suggest 
that the form of the diffusion equations remain the same 
for arbitrary momentum dependence. 

In particular the form of Eq. (3.25) proves the full 
analogy between singlet-triplet and charge-spin coupling 
in diffusive systems, [cf. Eqs. (2. 5-2. 6 )]. In Ref. [46], 
the analogy between the diffusion of spin in normal sys¬ 
tems and the triplet components was discussed. Here 
we can extend this result and find that the tensor C %, re¬ 
sponsible for the SHE in normal systems, is an additional 
source for the singlet-triplet conversion and, as we will see 
in the next sections, is at the root of magneto-electric ef¬ 
fects and the anomalous phase. Equations (3.22 -3.23) 
and (3.24-3.25) are the central equations of this work, 
which we now solve for different situations. In section 
VIB we go beyond this linear approximation. 










IV. THE EDELSTEIN EFFECT IN BULK 
SUPERCONDUCTORS FOR T -> T c 


In order to illustrate the usefulness of the SU(2) covari¬ 
ant quasiclassical equations presented above, we study 
here the magneto-electric effect and its inverse in bulk su¬ 
perconductors with an intrinsic SOC linear in momentum 
and derive the response coefficients in (1.3) and (1.4). 

We assume that the superconducting order parameter 
A is constant in magnitude but has a spatially dependent 
phase A (r) = |A|e 1( ^ r \ where X7p is assumed to be a 
constant vector. 

Let us first consider a diffusive superconductor. From 
(3.24) in the lowest order of S7p one obtains 


fs » 

k o n \ 


(4.1) 


and hence one can easily obtain the lowest order correc¬ 
tion to the triplet component from (3.25) : 


ft = 


IwJ 2m 


sgn (w„) 


(f + 4)' 


ab 


3 b jdj<p. (4.2) 


From Eq. (3.21) it becomes clear that the spin density 
is determined by the product of the singlet (4.1) and 
triplet (4.2) components which results in s a = Xidup 
with 


x“ = 4tt1Vo-^T Y, K 
2m 



(4.3) 


This is the Edelstein result generalized for arbitrary lin¬ 
ear in momentum SOC. 

With the help of Eqs. (3.24-3.25) we can also describe 
the inverse EE, the so-called spin-galvanic effect. We 
now assume a finite and spatially homogenous Mg and 
a zero phase gradient. In such a case one can obtain f t 
directly from Eq. (3.25), which is now proportional to 
Mg. By substitution of this result into the expression 
for the current, Eq. (3.20), and by noticing that only 
the second line contributes to the current we obtain ji = 
eygMg> with Xi given by Eq. (4.3) in agreement with 
Onsager reciprocity. 

In short, we are able to derive in a few lines the tensor 
(4.3), which describes the EE and inverse EE in super¬ 
conductors. Moreover, the expression (4.3) is valid for 
arbitrary linear in momentum spin-orbit effect and gen¬ 
eralizes the result obtained in Ref. [80] for the particular 
case of a Rashba SOC. If one assumes the same here, i.e. 
A y x = — a = —Ay and all the other components equal to 
zero, one obtains from Eq. (4.3) 




ay 


5%) 4ttNq 


Dt 

2m 


T E 

Wn> 0 


A 2 a 3 
ujI 2 \cu n \ + Da 2 


(4.4) 

that coincides with the expression obtained in Ref. [80]. 

If we neglect in Eq. (4.2) the Dyakonov-Perel relax¬ 
ation, then the triplet component is simply proportional 


to / t a ~ Aq (r). By substituting this into the expression 
for the current Eq. (3.20) one can easily show that 

<«> 

uj n > 0 n “ 


This expression suggests that a spatially inhomogenous 
magnetization together with SOC may also induce a fi¬ 
nite supercurrent. In this case the spin-galvanic effect 
scales with the square of the SOC parameter, in contrast 
to the a 3 dependency found previously for spatially uni¬ 
form magnetization. 

The same effects can be explored in the pure ballis¬ 
tic limit, for which Eqs. (3.22-3.23) apply. The singlet 
component in the lowest order in the SOC is given by 


fs 



1 -i 


vfth 

2uj n 



(4.6) 


whereas the triplet component can be obtained easily 
from Eq. (3.23) 


ft=~ 


A v F 

2 |cI (jjrp 


(v F rikPk + 2 c 


-l 


I ab 


rh? b ija 

- ~On (P. 

2m 3 


(4.7) 

By using Eq. (3.21) we obtain the Edelstein result 
S a = X^djP but now for an arbitrary linear in momen¬ 
tum SOC 


x a = _2,^rx 


2m 


E- 

w n >0 


A 2 


(v F n k Pk + 2 c 


-l 


I ab 




(4.8) 


Identically, we find ji = cx^Aq. 

In the particular case of a 2D systems with Rashba 
SOC we recover the Edelstein result for a ballistic super¬ 
conductor [79]: 


Xi = 


ttNqA 2 
4 v F m 


t y — 

X | 13 

u n > 0 \Wn\ 


(■ v F a ) 3 

(2cj n ) 2 + (v F a) 2 


(< 5 ; 


ay 



(4.9) 

The agreement between our and Edelstein results proves 
the validity of the expression (3.21) in the linearized ap¬ 
proximation. 

To conclude this section we note that for Rashba SOC 
in both cases, diffusive (4.4) and ballistic (4.9), Xi is 
proportional to a in the strong SOC limit (see also [134]), 
and to a 3 for weak spin-orbit interaction (see also [92]). 
So, the quasi-classic formalism is able to recover in an 
elegant way some well established results obtained after 
cumbersome diagrammatic [79, 80], and it also allows 
some easy generalizations of them. 


V. MAGNETO-ELECTRIC EFFECTS IN 
DIFFUSIVE JOSEPHSON JUNCTIONS 

We now turn to the central topic of the present work 
which is the description of magneto-electric effects in S- 
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X-S Josephson junctions and demonstrate their connec¬ 
tion to the anomalous phase problem. We first consider 
the diffusive limit and postpone the discussion of ballistic 
junctions for the next section. 

In particular we consider a S-X-S Josephson junction 
with an interlayer X of length L. We assume that the 
magnetic interactions are only finite in X and vanish in 
the S electrodes. Moreover, we assume that the struc¬ 
ture has infinite dimensions in the y — z plane and there¬ 
fore the GFs only depend on the x coordinate. The 
superconducting bulk solutions in the leads are written 
as /l = /Bcse _i ^ /2 and f R = /Bcse iv/2 , in the left 
(x < —L/2) and right (x > L/2) electrodes respectively, 
with 


/bcs 


+ A 2 


(5.1) 


whereas the normal metal fills the region —L/2 < x < 
L/2. 

We will consider both the highly resistive and the per¬ 
fectly transparent interfaces between the S and X parts. 
When the barrier transparency is low, the linearized ap¬ 
proximation is justified for all temperatures, whereas for 
transparent barriers, one is limited to temperatures close 
to the critical temperature of the junction. 

For the particular case of Rashba SOC in the X re¬ 
gion and an in-plane exchange field the Josephson cur¬ 
rent has been calculated in Ref. [92]. It has been 
shown explicitly that the current-phase relation is given 
by I = I c sm(ip — cpo). The anomalous phase ipo was cal¬ 
culated as a function of the strength of the spin fields, the 
temperature and the junction parameters. Here instead 
we focus in a generic linear-in-momentum SOC and we 
derive the expressions for the anomalous Josephson cur¬ 
rent in the lowest order of the spin fields. This will allow 
us to understand the link between the inverse EE and 
the (/?o-junctions. 


A. Diffusive junction with low transparency 
interfaces 


We first consider a S-X-S diffusive Josephson junction 
with highly resistive S-X interfaces. In this limit the lin¬ 
earization of the quasiclassical equations is justified for all 
temperatures. Our goal here is to determine the Joseph¬ 
son current through the junction, which in the linearized 
regime is given by Eq. (3.20). The components of the 
condensate function f s entering this expression, has to 
be obtained by solving the system (3.24)-(3.25) in the 
normal metal which couples the singlet with the triplet 
component. For the specific S-X-S geometry considered 
here this equations read: 


( 9 X - K l) fs+ sgn(w„) 


da - 


n = 0 Q ; 


(5.2) 


(vw// a - k 2 J?+ 

+ sgn(w„) (f ^ fs= 0 , (5.3) 

and the boundary conditions for the resistive interface 
(c/. Eq.(3.18)): 

(d x f s + sgn(w„) ^-d a jf) = ±7 f R>L 

\ 2m J x=±l /2 

dtf\x=±L/2 = 0 . (5.4) 

The expression Eq.(3.20), can be simplified by calcu¬ 
lating the current at the right interface (x = L/2) and 
by using the boundary condition (5.4): 

jx = ieirDNoT^ £ IV {fJ R - f s f R } x=L/2 ■ (5.5) 

^n >0 

It is clear from this equation that the correction to the 
current due to the spin-fields (the anomalous current) is 
proportional to Im [f^Sf s (L/2)\, where Sf s is the first 
correction to the singlet component due to the gauge po¬ 
tentials. In the absence of a phase difference between 
the S electrodes /# is real and the anomalous current is 
proportional to the imaginary part of the singlet compo¬ 
nent. According to Eq. (5.2), in the absence of spin-fields 
(exchange and SOC), there is no triplet component and 
the singlet component is real. Therefore no supercurrent 
flows at zero phase difference. 

In the presence of spin-fields there are two sources for 
singlet-triplet conversion, as seen from the second term 
in the l.h.s of Eq. (5.3). The first one is the exten¬ 
sively studied mechanism for singlet-triplet conversion in 
S/F junctions via the intrinsic exchange field Aq [1? ]. 
Inclusion of SOC leads to an additional singlet-triplet 
conversion mechanism described by the last term in the 
l.h.s of Eq (5.3). As discussed in section II, the singlet- 
triplet conversion in this case corresponds to the charge- 
spin conversion in normal systems with SOC. Conversely, 
once the triplet component is created, both mechanisms 
will convert it back to singlet, as can be seen in Eq.(5.2). 

The singlet-triplet-singlet conversion at the lowest or¬ 
ders in perturbation with respect to the spin-fields is 
schematized in Fig. 1. The black arrows represent the 
singlet-triplet conversion due to the exchange field which 
implies a tt /2 phase shift due to the i factor in front of Mo 
in Eqs. (5. 2-5. 3). The red arrows represent the singlet- 
triplet conversion due to the SOC, specifically due to the 
coupling term in Eqs. (5. 2-5. 3) proportional to No 

additional phase is associated with this latter process. If 
one follows the black path, i.e. the singlet-triplet-singlet 
conversion only due to the exchange field, the resulting 
contribution to the singlet component acquires a minus 
sign (a 7r shift) and it is proportional to Aq. This means 
that there is no anomalous phase 0 < ipo < n induced and 
hence no Josephson current flows when tp = 0. Similarly, 
if one follows the red path the resulting singlet compo¬ 
nent also remains real with no change of sign. From Fig.l 
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it becomes clear that a nontrivial <po only appears from 
the "cross-term” path that consist in one black and one 
red arrow. In other words, the mutual action of exchange 
field and SOC leads to a finite cpo and hence to a super¬ 
current even at zero phase difference. In this case the 
contribution to this current in the lowest order of the 
spin fields, is proportional to A^^difs between the ex¬ 
change field and the spin-current tensor, as anticipated 
in the introduction. 

In order to quantify this effect and calculate ipo in the 
S-X-S junctions it is convenient to introduce the singlet 
and triplet propagators associated with Eqs. (5. 2-5. 4): 

( d x - K l) K S (x, x') = -6(x- x') 

d x K s (x,x')\ x _ ±L/2 = 0 , (5.6) 


and 





k t (x, x') 


—S (x — x') 


(d x + P x )k t (x,x’) = 0 . (5.7) 


Thus, Eqs. (5. 2-5. 4) can be re-written as a set of integral 
equations 


fs (x) = / s (0) (x) - ^ [K s (x, L/2)f t a - K,{x t —L/2)} 


sgn (uj n ) / dxiK s (x,xi) x 

J -L/2 

ft( x i) , (5.8) 


i '^ Ks 



Figure 1: Schematic representation of the singlet-triplet- 

singlet conversion process at the lowest order with respect 
to the spin-fields. Black arrows represent the action of the 
exchange field, whereas red arrows encode the effect of the 
singlet-triplet coupling term due to the SOC. Only mixed red- 
black paths lead to the appearance of an anomalous phase 
<po in the singlet component and hence to a supercurrent in 
a S-X-S junction even without a phase bias between the S 
electrodes 


and 


ft ( X ) = Sgn 



dx i x 


Kf (x, Xl ) 



+ 2^- 


fs {xi) 


(5.9) 


f L J and 

the second term in Eq. (5.8) takes into account the 
boundary condition (5.4) 

The K s propagator can be obtained from Eq. (5.6) 


Here /j 0) = 7 ( K s ( *, | ) f R + K s 




K s (xi,x 2 ) = 

coshK^ (L - \x\ - x 2 \) + coshK^ (xi + x 2 ) , , 

2kcjS inh/q^L 

whereas the equations for the triplet kernel, Eqs. (5.7), 
can be written in the form of an integral equation which 
is convenient for the subsequent perturbative analysis: 


K t (x 1 ,x 2 ) = e PxXl K s (xi,x 2 )e PxX2 + 

r L / 2 - . „ 

— e~ PxXl / dyK s ( Xl ,y)e Px yr ± K t (y,x 2 ) (5.11) 

J-L/2 

where T± = — Py — P%. 

In the lowest order of the gauge potentials one can 
obtain the correction Sf s to the singlet component by 
substituting the result (5.10) into Eqs. (5.8)-(5.9). We 
consider here only the "cross-term” correction Sf s pro¬ 
portional to both the exchange field Aq and the spin- 
current 3i and which is responsible for the anomalous 
phase-shift: 


Sf s (L/2) = iA a 0 d b x A^x 


2mD 


/ L/2 /.L/2 

dyi / 

-L/2 J- 


—L/2 


dy 2 Kf(y 2 , yi ) 


cosh [k u (yi -y 2 )} 


, sinh KujL 


(5.12) 


In principle, one has all the elements to solve Eqs. (5.8- 
5.9), for example recursively by performing a perturba¬ 
tive expansion in the gauge potentials. Here, in order to 
get analytical compact expressions we restrict our analy¬ 
sis to the short junction limit, i.e. L min(^J 1 , \ Ak\~ 1 )- 
In this case K s « /qJ 2 L -1 ( cf\ Eq. (5.10)) and from Eq. 
(5.11) it is easy to verify that K t reads 


Kt 


(4 + F) 

L 


(5.13) 


We are interested in calculating the anomalous phase 
(po which can be obtained by noticing that the current 
(5.5) can be written as 


jx = jc sin (<p - <po) W jc sin <p - (p 0 j c cos tp (5.14) 












for a small Lpo. The anomalous phase ipo can be obtained 
by setting ip = 0 and dividing by the critical current j c 
in the absence of SOC. In the short junction limit this is 
given by: 


jc = AenDNoTc'r 2 £ (5.15) 

' hvf. 1 j 

W n >0 W 

We follow this procedure and from Eq.(5.5) and Eqs. 
(5.12-5.13) we obtain 


X 


/bcs 


^0 


2 mD 


L- 


x>0 


na 

^0 


(*l + r±)' 


UJ. 


E : 

>0 


BCS 


ab 


(5.16) 

This expression clearly shows the relation between the 
appearance of the anomalous phase, cpo, and the inverse 
Edelstein effect in bulk systems. Both, the Josephson 
current (proportional in the linearized case to tpo) and the 
bulk supercurrent are proportional to Ao8 x , ^.e. both are 
generated from the mutual action of the exchange field 
and the SOC. 

It is worth noticing that in the present case of low 
transparent interfaces, the anomalous phase grows lin¬ 
early with L, the length of the junction (5.16). In the 
next subsection we show that in the case of a transpar¬ 
ent barrier the anomalous phase behaves like L 3 . 

In the particular case of a 2D situation, with a SOC 
coupling of Rashba (described by the parameter a) and 
Dresselhaus (/?) type we obtain from Eq. (5.16): 


/bcs ~ a ^o) (^ 2 ~ fi 2 ) 
tl£z 0 < 2uJn + D(a 2 ^P 2 ) 

2m \ ' /bcs 

2 -^ ac.2 


(5.17) 


Besides controlling the anomalous phase and hence the 
Josephson current by tuning the external magnetic field, 
this expression also suggests that the current can be con¬ 
trolled by tuning the Rashba SOC by means of an exter¬ 
nal gate. In the particular case that a = (3 the anomalous 
phase is zero and no supercurrent will flow. 


B. Diffusive junction with transparent interfaces 


We now briefly consider the limit of a full transpar¬ 
ent barrier. In that case one assumes continuity of the 
quasiclassical GFs at the S-X interfaces. The problem 
is then formally the same as in the previous section, ex¬ 
cept that the second equations in (5.6) and (5.7), for the 
propagators K s and K t are replaced by: 


Kt (^>^2) 


= 0 

xi=±Z//2 


(5.18) 
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respectively. In this case one should re¬ 
move the second term in Eq. (5.8) and 

/i 0) (x) = /z, sinh (L/2 — x) / sinh (/q^L) + 

f R sinh (L/2 + x) / sinh (k u L). 

Now the singlet propagator is given by: 


K s (xi,x 2 ) = 

cosh Km (xi + X 2 ) - cosh Kv (L - \xi - X 2 \) .g 
sinh 

In the short junction limit K s is proportional to L and 
it is temperature independent. From Eq. (5.11) K t ~ 
K s . Thus, in this case the anomalous phase-shift is also 
temperature independent and proportional to 


T T 3 


(5.20) 


In contrast to the case of finite barrier resistance, Eq. 
(5.16), the anomalous phase scale with L 3 . This means 
that in short junctions a finite barrier resistance be¬ 
tween the S and the normal metal favors the growth of 
ipo. These results generalize those presented recently in 
Ref. [92] for the particular case of Rashba SOC. 

We can then conclude that the anomalous phase, at 
lowest order in the gauge potentials, is proportional to 
Aq8x: independently of the type of interface. 


VI. MAGNETO-ELECTRIC EFFECTS IN 
BALLISTIC JOSEPHSON JUNCTIONS 

In this section we consider a pure ballistic S-X-S junc¬ 
tion,^ e. we solve (3.8) in the limit r —» 00. As before, the 
junction is along the x-axis and the two superconducting 
electrodes at position x < —L/2 and x > L/2. The spin 
fields, both exchange and SOC, are only finite in the X 
region. We assume that the the transverse dimensions 
of the junction are very large, and therefore the GFs de¬ 
pends on xand only weakly on y, z. We also assume that 
the interfaces between X and S are perfectly transparent. 

In the next subsection we first analyze the Josephson 
current for temperatures close to the superconducting 
critical temperature T c , and make a connection with the 
diffusive structures studied in the previous section. In 
the second subsection we derive analytical expressions 
for the anomalous current at arbitrary temperature for 
the case of small spin fields. 


A. Ballistic junction at T —»• T c 

In the case of large enough temperatures we analyze 
the linearized Eilenberger equation The solutions for 
the singlet and triplet components in equations (3.22) 
and (3.23) can be written as propagation in two di- 
rections f s ,t (~L/2 < x < L/2) = f> t (x) 0 (uj n /n x ) + 
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ff,t ( x ) e (-u n /n x ) with 


B. Arbitrary temperatures 


< _ A (Lf 2) 2uj n (x-L/2)/v F n x , f 

h M Jl 


dy 


L/2 V F^ X 

e _ 2oJ n (x—y),VFn x ( Ug + iXI JL) /»< (y) (6.1) 


and 


rP = f 

Jl 


dy 


2a } n {x—y)/vF'n x 


L /2 v F^x 
PiUi (x—y)/n 


f 


( 6 - 2 ) 


2m drid 


In the opposite propagation direction ff t are found from 
ff t by substituting L/2 —L/2. 

In analogy with the diffusive case (c/. Fig 1), expres¬ 
sions (6.1) and (6.2) show explicitly the effect of the SOC 
on the condensate function. In the absence of SOC the 
exchange field Ao is the only source for singlet-triplet 
conversion. The manifestation of the triplet component 
in S-F-S junctions has been extensively studied in the 
past (see [1, 2] for reviews). As discussed in section II, 
the imaginary unit i in front of the Ao terms leads to a 
7r/2 phase shift. In the case of a finite SOC the gauge- 
field, is an additional source of triplet correlations. 
Notice that in the ballistic case, $ij not only couples 
the singlet and triplet components, but also the s-p- wave 
components of the condensate [143]. Moreover, the term 
e - p i n ix/n x (6.2) leads to a momentum dependent ro¬ 
tation of the triplet component in the spin-space. 

The origin of the anomalous phase cpo can be easily un¬ 
derstood in the lowest order in the spin fields. Assuming 
a vanishing phase difference between the superconduc¬ 
tors and combining Eqs.(6.1)-(6.2) with the expression 
for the current (3.19), one obtains for the first nontrivial 

( \ ab 
e -Piriix/n x \ j^b ^ 

'?ijPj' b A' 0 - This correction is proportional to S^Aq and 
coincides with those obtained in bulk superconductors 
with SOC (section IV) and a diffusive S-X-S junctions 
(section V). 

Quantitatively, a compact analytical solution for the 
current at zero-phase difference can be obtained from 
Eqs. (6. 1-6.2) in the short junction limit, i.e. for 
L <C Vp/^UOn'. 


The previous result for the current has been obtained 
at temperatures close to the critical one. We now con¬ 
sider an arbitrary temperature and calculate the current 
up to the lowest order in the gauge field 9^. In order 
to calculate the current from Eq.(3.10) to the lowest or¬ 
der in 3^ij we need to compute the first two components 
matrix g = g^ + + • • • . 

At zeroth order in 5T i<7 - the ballistic equation reduces to 


v F n x 


dg 

dx 


i r 3 (icj n + M 0 ) + v F rijAj + A,g^ 


which admits for solution 


(6.4) 


g(°) = u(x)g^ 0) u 1 (x)+g OQ: (6.5) 


with g^ a constant matrix found from the boundary 
conditions. The propagator u (x) in Eq. (6.5) is given by 


u (x) = exp 


. t 3 {icd n + Ao) + vpfijAj + A 

v F n x ^ 


( 6 . 6 ) 


when we assume that neither Ao nor A nor Aj depend on 
the position, ^describes how the function go “propagates” 
from its value at x = 0, ? t° any point x. 

The constant g^ in (6.5) satisfies 

[r 3 (icj n + Ao) + v F rijAj + A, g^ \ = 0 , (6.7) 

and describes the bulk contribution deep inside the su¬ 
perconductor. Notice that according to Eqs.(6.7)-(6.6), 
[<?oo, v] = 0, hence g^ can not be obtained by the appli¬ 
cation of (6.6). As we will see below [see (6.13)], the so¬ 
lutions of u in a superconductor are evanescent waves, so 
the contribution ug^vT 1 vanishes deep inside the super¬ 
conductor, whereas the contribution g^ remains finite. 

The first-order correction with respect to the gauge- 
field satisfies 


v F n x 


dgW 

dx 


niPii dg (0) \ 

2m ’ drij ) 


r 3 (iw n + A 0 ) + vpnjAj + A, g (0) 


( 6 . 8 ) 


and so 


3x (¥> = 0) = _*/a 2 T c J^«x 


E f 


X 

UJ n > 0 


0 -2u n L/v F \n x \ 


1 + 2^- 


. (6.3) 


where = d!JMg. Thus the anomalous current 

is generated by the spin-polarization Mg via the spin- 
current This is the spin-galvanic effect, discussed in 
the previous sections, for a ballistic S-X- S junction. 


g (1) = U (x) 5o 1} (■ x ) u 1 (x) , (6.9) 

with a position-dependent p 1 matrix, which reads 

5o 1} = 9l 


f dz , t f ntfij dgW \ . 

+ / - « W 

Jo v Fn x { 2m drij ) 


v F n x 

where gi is a constant matrix. 


, ( 6 . 10 ) 
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The current can be written in powers of 3^, j x = 
jl 0) +j* 1) H-with [see (3.10)] 


.( 0 ) _ ienN 0 v F 


X Tr 


*(o) 


T3 


( 6 . 11 ) 


and the first order correction 




( 6 . 12 ) 


Notice that the second line in (6.10) vanishes after the 
angular average. We then need to obtain g^ and g\ 
to determine the current through the S-X-S Josephson 
junction. 

We separate the solution of the problem in the three 
regions: the two superconductors (x > L/ 2 and x < 
—L/2) and the normal region (-L/ 2 < x < L/ 2 ). One 
can check that in the superconductors: 


H x <~\ 


L 


9 ( x > - ) = e 


* \s L g L T+S{ 


+ 9oo 


e 1T3 4 , 


r_5L+5oo e 1T3 % 


with r± = (ri ± ir 2 ) /2 and 

TsUJn +T 2 a r 3 sinh r) + r 2 


3oO — 


+ A 2 


cosh 77 


(6.13) 


(6.14) 


where f2o = u (A = 0) [see ( 6 . 6 )] is a unitary matrix that 
can be written as 


uo (x) = e ~ UJriT3X / VFnx 


u 0 \ 

0 ) ' 


The spin matrices u and u are defined as 


u (x, n) = exp 


. A 0 T vpifinAj 

1 -- —-x 


v F n x 


u (x, n) = cr y u * (x, — n) 


= exp 


. — Aq + vprijAj 


v F n x 


(6.18) 


The matrices go and g\ in Eq. (6.16) are obtained from 
the boundary conditions, assuming continuity of the GFs 
at the left and right boundaries. At zeroth order we 
obtain 


9o 


9 o_ fo A 

-/o —go J 


(6.19) 


go 



2 cosh 2 y + Tr {U (L)} 


+ 2 sinh 2 y 


( 6 . 20 ) 


e x g (|) +e -x g (-|) 
J0 2 cosh 2 y + Tr {U (L)} 


( 6 . 21 ) 


Sl,r 


e^l + irie 11 ! ^ T3 CO sh?7(cc±L/2)/^o 
a/2 cosh 77 


(6.15) 


X = —--I- arcsinh ^ + i— (6.22) 

A 2 


where sinh 77 = huj n /A , and £0 = hv F /A is the su¬ 
perconducting coherence length. The matrices gL,R ~ 

+ • • • have been expanded in power of 3^- ; 

# 1 °# ^ are constant matrices found from boundary con¬ 

ditions order by order, g^ is present at the zeroth order 
only. 

In the normal region, the solution reads 

9 < X < = Mo (x) g 0 ul (x) + 

+ Mo (x) 9o 1} ( X) ul(x) + ■■■ , (6.16) 


U (x) = u (x) u (—x) (6.23) 

whereas U (x) = u (x) u {—x) is its time-reversal conju¬ 
gate. We here give only the contribution corresponding 
to the positive projection of the Fermi velocity, the neg¬ 
ative projection can be found straightforwardly. 

The matrices entering the first-order correction, Eq. 
(6.17), have the following form in Nambu space 

* = ( -h 4 ) a,,tl ® “ ( 4 -e ) ' < 6 - 24 > 

with 


where 

5o 1} (•■'■) </i + / © (z) dz , (6.17) 

Jo 


&(z) 


1 

v F n x 


4 (z) 


nSFjj dup (z) gpUp (z) 
2m ’ drij 


Uo{z) , 



\ r L / 2 ^ r L /2 

+ 7; © • 9odz - - S • f 0 dz . (6.25) 

Z J-L/2 1 J — L/2 

After multiplication by n x and taking the angular average 
the first line of this equation vanishes. The second line of 
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(6.25) can be simplified using the normalization condition 
gl - /o/o = 1 available for the zeroth order correction. 
We obtain 


= ifI N£vr £ f 1/2 dz £ 

'• ^nJ-L/2 a =± 


Tr 


2 

ntfn 


Wn ^>0 


'f„3£-(„|t| 8g n(„,)) , (6.26) 


2m 


* dn« 


" On, 


with 


f± (*) = -2i- 


e ±x± [7(2;+—)+ e Tx± U \ z — — 


2 cosh 2y± + Tr {1/ (L)} 


CoU L Ur, (0 , 

X± = -i—I + arcsinh — ± (6.27) 

and f± (z, n) = cr y f± (z, — n) a y its time reversal conju¬ 
gate. 

If Aq commutes with Ai, then u(x)u(—x) = 

exp [iAqx/v f^x] becomes independent of the SOC [see 
the definitions (6.18)], and the contribution (6.26) van¬ 
ishes. Therefore we expect (6.26) to be proportional to 
$ij [AojAj] oc TijToj at the smallest order in the gauge- 
fields. 

By expanding the expression (6.26) in the gauge- 
potentials, up to the term proportional to the electric-like 
field one obtains 



Figure 2: The temperature dependence of the averaged 
anomalous phase shift (6.33) in a log-scale In (^o (T) /4>o), 
where 4>o = — hL 3 Tr {T xy T y o} /6 Ef- We have approximated 

A (T) « 1, 764T C tanh ^1, 74y / T c /T — l^j which is the usual 
interpolation for the temperature dependence of the super¬ 
conducting gap. The curves are given for different ratio 
of L/fo = {0.01, 0.06, 0.11, 0.16}, with £ 0 = Hvf/ A 0 and 
Ao = 1.764T C the gap at zero temperature. Note that ipo 
does not vanish when T —)> T c . 


f ±(z) 


—(A) 

2 cosh 2 x± \ v x \ 2 / 

/ 2L4q 

2 cosh 2 x± \ Vx 



(.z + L/2) A 


Vk [Ak, A] 


L\\ 2 (z-L/2) 2 


Vk [Ak, A)] 


(6.28) 


where Tr {U (L)} « 2 in the small gauge-field limit. Be¬ 
sides the terms proportional to Ao-only, responsible for 
the oscillations of S/F proximity effect, the SOC Aj only 
appears in the electric-field construction (the last term 
on each line), due to symmetry with respect to the time- 
reversal. After angular averaging only the last contribu¬ 
tions of the two lines are non-zero. This leads to 



jo 


-2 


L 3 

3 Ep 


x 


TY{AiAo} 



n j 

1 + 2 —|- 
n± 



(6.30) 


J —=4(\n x \M) , 
Jo 


(6.29) 


with jo = tt evpNoT and 

M = 9 / tanh ( l PT. 4 . arcsinh + i— 

jfeo l \v F \n x \ A 2 

(6.31) 

(note that the sum over j applies inside the angular aver¬ 
aging as well). As in all previous examples the anomalous 
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current is proportional to Tr {J'xk^ko} = 9 x Aq, where 
the later form suggests our expressions are valid beyond 
the linear-in-momentum-SOC approximation, given any 
spin current 3f and paramagnetic interaction Aq. 

Close to the critical temperature M « 
A 2 sin(/9^ u; >0 e~ 2ujL / VF \ n X \/2u 2 and we recover 
( 6 . 3) . 

Commonly, the concept of a po ~]unction is defined for 
junctions with a sinusoidal current-phase relation. This 
is valid at temperatures close to the critical temperature 
or in the case of a weak proximity effect between the S 
electrodes and the X bridge. However, in several cases 
the current-phase relation is more complex and higher 
harmonics are involved [132]. This is the case of the bal¬ 
listic junction studied here with a current-phase relation 
given by sum of Eqs. (6.29) and (6.30). In such cases the 
Po is defined as the phase difference across the junction 
that minimize the energy, or equivalently, as the phase 
difference imposed to the junction in order to get a zero 
current state, i.e. j (po) = 0. In our perturbative analy¬ 
sis ipo is small and hence 


<Po = 


l (y = o) 
M<„=o ' 


(6.32) 


It is clear that j (ip = 0) = j x \ whereas cl P j _ 0 = 


d v j. 


; ( 0 ) 


c p=0 


and from Eqs. (6.29) and (6.30) we obtain 


<P o = 


hL 3 


n-an-a A /it 

difi ' ' 



ip=0 


6Ep 




ip =0 


(6.33) 

In Fig. 2 we show the temperature dependence of 
Po for the ballistic junction for a 2D system when only 
^xy^yo non-zero. We assume a circular Fermi surface, 
n x = cos 6 and n y = sin#. We plot the anomalous phase 
for different junction lengths. 


VII. DISCUSSION AND CONCLUSIONS 


In order to verify our findings and prove the existence 
of the anomalous po phase one can design a supercon¬ 
ducting ring interrupted by a semiconducting link with a 
strong SOC, similar to the one used recently in Ref. [144] 
for the characterization of the current phase relation of a 
Nb/3D-HgTe/Nb junction or in [145] for the observation 
of a spontaneous supercurrent induced by a ferromag¬ 
netic 7r-junction. A schematically view of the proposed 
setup is shown in Fig. 3: It consists in a superconducting 
ring (green) grown on top of a semiconductor or a metal¬ 
lic substrate with strong SOC (grey). In order to isolate 
electrically the S ring from the semiconductor one can for 
example add an insulating barrier (blue) under the ring. 


If a magnetic field is applied in the plane of the ring, 
it will act as a Zeeman field and hence, according to our 
previous results, it will create a spontaneous circulating 
supercurrent, see [7, 146] for more details. This super¬ 
current will generate a magnetic flux that in principle 
can be measured by a second loop [144] or a micro-Hall 
sensor [145]. 

In the case when the bridge is made of a 2D semicon¬ 
ductor with a generic SOC described by a combination of 
Rashba and Dresselhaus terms: A x = —aa v + /3a x and 
A y = acr x — /3a y the generated supercurrent should be 
proportional to 


js oc (a 2 - P 2 ) (h x p + h y a) . 

Thus the current depends on the direction of the applied 
magnetic field. In particular, for a field perpendicular 
to the 2D gas the effect should vanishes. In addition 
by applying a gate voltage one could modify the ratio 
between Dresselhaus and Rashba interactions and hence 
control the supercurrent flow. We thus expect that the 
dependency of the spontaneous supercurrent with respect 
to the orientation of the magnetic field and/or the gate 
voltage realizes a clear demonstration of the spin-galvanic 
effect in Josephson systems. 

Instead of using a semiconducting bridge one could 
grow the superconducting loop on top of a metallic sub¬ 
strate. Metals with strong SOC, like Pt and Ta, are good 
candidates to observe the po~j unction behavior, but also 
an ultra-thin layer of Pb might be used [147]. In such 
a case probably one cannot control the (/?o-shift using a 
gate, but a spontaneous circulating current might still be 
controlled by switching the in-plane external field on and 
off. 

Eventually, the existence of a magneto-electric phase- 
shift po can be probed by measuring the Shapiro steps 
in S-X-S Josephson junctions as suggested in Ref. [148]. 



Figure 3: One possible experimental setup to verify the gen¬ 
eration of a spontaneous current in S/X/S Josephson junction. 
It consists of a truncated superconducting loop (in green) de¬ 
posited on top of a material exhibiting SOC (in grey). In 
order to isolate electrically part of the S loop from the con¬ 
ducting substrate we assume an insulating layer between them 
(in blue). By applying an in-plane magnetic field a circulat¬ 
ing supercurrent might be generated which in turn induces a 
magnetic flux that can be measured via an extra pick-up coil 
(shaded grey). For more details see discussion in the main 
text. 
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In conclusion, we have demonstrated that the inverse 
Edelstein effect, also called spin-galvanic effect, and the 
appearance of an anomalous phase-shift cpo in Joseph- 
son junctions are the two sides of the same coin. We 
presented a full SU(2) covariant quasi-classic formalism 
that allows to study these magneto-electric phenomena 
in bulk and hybrid superconducting structures with ar¬ 
bitrary linear-in-momentum SOC (section III). 

With the help of our quasi-classic transport formal¬ 
ism we derived the Edelstein effect close to the critical 
temperature of a bulk superconductor, recovering the 
Edelstein’s result in a very compact way (section IV) 
and generalizing it for the case of an arbitrary linear- 
in-momentum SOC. We have shown that the Edelstein 
effect and its inverse are reciprocal in the sense of the 
Onsager relations, both in ballistic and diffusive super¬ 
conducting systems: A static supercurrent can induce a 
finite magnetization due to the presence of a spin-orbit 
coupling, and reciprocally a finite magnetization pro¬ 
duces a finite supercurrent in a bulk system. We have 
demonstrated that the linear-response tensor is directly 
proportional to the equilibrium spin-current tensor d ?. 

We have also generalized this result to inhomogeneous 
systems. In particular we have studied the current-phase 
relation of a Josephson junction consisting of two super¬ 
conductors coupled via a normal metal with both SOC 
and spin-splitting field. We have demonstrated that a 
supercurrent can flow even if the phase difference be¬ 
tween the S electrodes is zero. This current is associ¬ 
ated to an anomalous phase-shift This result holds 
for both ballistic (section VI) and diffusive systems (sec¬ 


tion V), for arbitrary linear-in-momentum spin-orbit cou¬ 
pling, and for arbitrary barrier resistance between the 
superconductor and the normal metal. For all these sit¬ 
uations we have demonstrated that SU(2) gauge-fields 
are the only objects of relevance in the phenomenology 
of the (^o-shift, and in particular we have shown that 
c po (x A^i = Aq j'fji’ ^* e * the anomalous phase-shift is 
proportional to the SU(2) electric and magnetic fields, 
or equivalently to the spin-current tensor. We thus di¬ 
rectly linked the anomalous phase-shift in superconduct¬ 
ing systems to the inverse Edelstein effect (also known 
as the spin-galvanic effect) extensively studied in normal 
systems. 


Acknowledgments 

F.K. thanks F. Hassler and G. Viola for daily stimulat¬ 
ing discussions, as well as A.I. Buzdin for insightful dis¬ 
cussions. F.K. is grateful for support from the Alexander 
von Humboldt foundation during his stay in the JARA- 
Institute for Quantum Information at RWTH-Aachen, 
Germany, where part of this work was done. The work of 
F.S.B. was supported by Spanish Ministerio de Econorma 
y Competitividad (MINECO) through the Project No. 
FIS2011-28851- C02-02 and the Basque Government un¬ 
der UPV/EHU Project No. IT-756-13. I.V.T. acknowl¬ 
edges support from the Spanish Grant FIS2013-46159- 
C3-1-P, and from the “Grupos Consolidados UPV/EHU 
del Gobierno Vasco” (Gant No. IT578-13) 


[1] A. I. Buzdin, Reviews of Modern Physics 77, 935 (2005). 

[2] F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Reviews 
of Modern Physics 77, 1321 (2005). 

[3] P. SanGiorgio, S. Reymond, M. R. Beasley, J. H. Kwon, 
and K. Char, Physical Review Letters 100, 237002 
(2008). 

[4] K. M. Boden, W. P. Pratt, and N. O. Birge, Physical 
Review B 84, 020510 (2011), arXiv: 1108.5613 . 

[5] T. Kontos, M. Aprili, J. Lesueur, and X. Grison, Phys¬ 
ical Review Letters 86, 304 (2001). 

[6] A. Cottet, Physical Review B 76, 224505 (2007), 
arXiv:0704.3975 . 

[7] L. Bulaevskii, V. V. Kuzii, and A. A. Sobyanin, JETP 
Letters 25, 290 (1977). 

[8] A. I. Buzdin, L. Bulaevskii, and S. V. Panyukov, Sov. 
Phys. JETP 35, 178 (1982). 

[9] V. V. Ryazanov, V. Oboznov, A. Rusanov, A. Vereten¬ 
nikov, A. Golubov, and J. Aarts, Physical Review Let¬ 
ters 86, 2427 (2001). 

[10] T. Kontos, M. Aprili, J. Lesueur, F. Genet, B. Stephani- 
dis, and R. Boursier, Physical Review Letters 89, 
137007 (2002). 

[11] F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Physical 
Review B 64, 134506 (2001). 

[12] F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Physical 
Review Letters 86, 3140 (2001). 


[13] J. W. A. Robinson, J. D. S. Witt, and M. G. Blamire, 
Science (New York, N.Y.) 329, 59 (2010). 

[14] T. S. Khaire, M. A. Khasawneh, W. P. Pratt, and 
N. O. Birge, Physical Review Letters 104, 4 (2009), 
arXiv:0912.0205 . 

[15] M. Anwar, F. Czeschka, M. Hesselberth, M. Porcu, and 
J. Aarts, Physical Review B 82, 100501 (2010). 

[16] M. S. Anwar, M. Veldhorst, A. Brinkman, and 
J. Aarts, Applied Physics Letters 100, 052602 (2012), 
arXivil 111.5809 . 

[17] I. T. M. Usman, K. A. Yates, J. D. Moore, K. Morri¬ 
son, V. K. Pecharsky, K. A. Gschneidner, T. Verhagen, 
J. Aarts, V. I. Zverev, J. W. A. Robinson, J. D. S. Witt, 
M. G. Blamire, and L. F. Cohen, Physical Review B 
83, 144518 (2011). 

[18] C. Klose, T. S. Khaire, Y. Wang, W. P. Pratt, N. O. 
Birge, B. J. McMorran, T. P. Ginley, J. A. Borchers, 
B. J. Kirby, B. B. Maranville, and J. Unguris, Physical 
Review Letters 108, 127002 (2012). 

[19] E. C. Gingrich, P. Quarterman, Y. Wang, R. Loloee, 
W. P. Pratt, and N. O. Birge, Physical Review B 86, 
224506 (2012), arXiv:1208.3118 . 

[20] J. D. S. Witt, J. W. A. Robinson, and M. G. Blamire, 
Physical Review B 85, 184526 (2012). 

[21] J. W. A. Robinson, F. Chiodi, M. Egilmez, G. B. Halasz, 
and M. G. Blamire, Scientific Reports 2 , 699 (2012). 



17 


[22] Y. Kalcheim, O. Millo, M. Egilmez, J. W. A. Robin¬ 
son, and M. G. Blamire, Physical Review B 85, 104504 
( 2012 ). 

[23] Y. Kalcheim, I. Felner, O. Millo, T. Kirzhner, G. Ko- 
ren, A. Di Bernardo, M. Egilmez, M. G. Blamire, and 
J. W. A. Robinson, Physical Review B 89, 180506 
(2014). 

[24] M. Alidoust, K. Halterman, and O. T. Vails, Physical 
Review B 92, 014508 (2015), arXiv: 1506.05469 . 

[25] N. Banerjee, J. Robinson, and M. G. Blamire, Nature 
Communications 5, 4771 (2014). 

[26] P. V. Leksin, N. N. Garif’yctnov, I. a. Garifullin, Y. V. 
Fominov, J. Schumann, Y. Krupskaya, V. Kataev, O. G. 
Schmidt, and B. Buchner, Physical Review Letters 109, 
057005 (2012), arXiv: 1207.0727 . 

[27] N. Banerjee, C. B. Smiet, R. G. J. Smits, A. Ozaeta, 
F. S. Bergeret, M. G. Blamire, and J. W. A. Robinson, 
Nature Communications 5, 3048 (2014). 

[28] X. L. Wang, A. Di Bernardo, N. Banerjee, A. Wells, 
F. S. Bergeret, M. G. Blamire, and J. W. A. Robinson, 
Physical Review B 89, 140508 (2014). 

[29] F. Chiodi, J. D. S. Witt, R. G. J. Smits, L. Qu, G. B. 
Halasz, C.-T. Wu, O. T. Vails, K. Halterman, J. W. a. 
Robinson, and M. G. Blamire, EPL (Europhysics Let¬ 
ters) 101 , 37002 (2013), arXiv: 1211.1169 . 

[30] Y. Gu, J. W. A. Robinson, M. Bianchetti, N. A. Stel- 
mashenko, D. Astill, F. M. Grosche, J. L. MacManus- 
Driscoll, and M. G. Blamire, Applied Physics Letters 
Materials 2 , 046103 (2014). 

[31] A. A. Jara, C. Safranski, I. N. Krivorotov, C.-T. Wu, 
A. N. Malmi-Kakkada, O. T. Vails, and K. Halterman, 
Physical Review B 89, 184502 (2014), arXiv: 1404.2304 

[32] A. Rusanov, J. Aarts, M. Hesselberth, and A. I. Buzdin, 
Physical Review Letters 93, 57002 (2004). 

[33] A. Rusanov, S. Habraken, and J. Aarts, Physical Re¬ 
view B 73, 60505 (2006). 

[34] M. Eschrig, Physics Today 64, 43 (2011). 

[35] J. Linder and J. W. A. Robinson, Nature Physics 11, 
307 (2015). 

[36] F. S. Bergeret, A. Verso, and A. F. Volkov, Physical 
Review B 86, 1 (2012), arXiv: 1211.0870 . 

[37] A. Pal, Z. Barber, J. Robinson, and M. Blamire, Nature 
Communications 5, 3340 (2014). 

[38] D. Massarotti, A. Pal, G. Rotoli, L. Longobardi, M. G. 
Blamire, and F. Tafuri, Nature Communications 6, 
7376 (2015). 

[39] F. S. Bergeret and F. Giazotto, Physical Review B 88, 
014515 (2013), arXiv: 1305.6301 . 

[40] F. S. Bergeret and F. Giazotto, Physical Review B 89, 
054505 (2014), arXiv: 1312.1945 . 

[41] A. Ozaeta, P. Virtanen, F. S. Bergeret, and T. T. 
Heikkila, Physical Review Letters 112 , 057001 (2014). 

[42] F. Giazotto, J. W. A. Robinson, J. S. Moodera, and 
F. S. Bergeret, Applied Physics Letters 105, 062602 
(2014), arXiv:1405.0565 . 

[43] F. Giazotto, T. Heikkila, and F. Bergeret, Physical Re¬ 
view Letters 114, 067001 (2015), arXiv: 1403.1231 . 

[44] U. Nasti, L. Parlato, M. Ejrnaes, R. Cristiano, T. Taino, 
H. Myoren, R. Sobolewski, and G. Pepe, Physical Re¬ 
view B 92, 014501 (2015). 

[45] F. S. Bergeret and I. V. Tokatly, Physical Review Let¬ 
ters 110, 117003 (2013), arXiv:1211.3084 . 

[46] F. S. Bergeret and I. V. Tokatly, Physical Review B 89, 


134517 (2014), arXiv:1402.1025 . 

[47] S. H. Jacobsen, J. A. Ouassou, and J. Linder, (2015), 
arXiv: 1503.06835 . 

[48] J. Arjoranta and T. T. Heikkila, (2015), 
arXiv: 1507.02320 . 

[49] M. Alidoust and K. Halterman, New Journal of Physics 
17, 033001 (2015), arXiv: 1504.05950 . 

[50] S. Maekawa, S. O. Valenzuela, E. Saitoh, and 
T. Kimura, Spin current (Oxford University Press, 
2012 ). 

[51] T. Kuschel and G. Reiss, Nature Nanotechnology 10, 
22 (2014). 

[52] A. Manchon, Nature Physics 10, 340 (2014). 

[53] C. Ciccarelli, K. M. D. Hals, A. Irvine, V. Novak, 
Y. Tserkovnyak, H. Kurebayashi, A. Brataas, and 
A. Ferguson, Nature Nanotechnology 10, 50 (2014), 
arXiv: 1411.2779 . 

[54] M. Dyakonov and V. Perel, JETP Letters 13, 467 
(1971). 

[55] M. Dyakonov and V. Perel, Physics Letters A 35, 459 
(1971). 

[56] J.-N. Chazalviel, Physical Review B 11, 3918 (1975). 

[57] J. E. Hirsch, Physical Review Letters 83, 1834 (1999), 
arXiv:cond-mat/9906160 . 

[58] M. Dyakonov and A. Khaetskii, in Spin Physics in Semi¬ 
conductors , Springer Series in Solid-State Sciences, Vol. 
157, edited by M. Dyakonov (Springer, Berlin, Heidel¬ 
berg, 2006) Chap. 8, pp. 211-243. 

[59] E. G. Mishchenko, A. V. Shytov, and B. I. Halperin, 
Physical Review Letters 93, 226602 (2004), arXiv:cond- 
mat/0406730 . 

[60] Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. 
Awschalom, Physical Review Letters 93, 176601 (2004), 
arXiv:cond-mat/0403407 . 

[61] Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. 
Awschalom, Science (New York, N.Y.) 306, 1910 (2004). 

[62] J. Wunderlich, B. Kaestner, J. Sinova, and T. Jung- 
wirth, Physical Review Letters 94, 047204 (2005), 
arXiv:cond-mat/0410295 . 

[63] R. Raimondi, C. Gorini, P. Schwab, and M. Dzierzawa, 
Physical Review B 74, 035340 (2006), arXiv:cond- 
mat/0601525 . 

[64] R. Raimondi, P. Schwab, C. Gorini, and G. Vignale, 
Annalen der Physik 524, 153 (2012), arXiv: 1110.5279 . 

[65] S. O. Valenzuela and M. Tinkham, Nature 442, 176 
(2006), arXiv:cond-mat/0605423 . 

[66] M. Morota, Y. Niimi, K. Ohnishi, D. H. Wei, T. Tanaka, 
H. Kontani, T. Kimura, and Y. Otani, Physical Review 
B 83, 174405 (2011), arXiv:1008.0158 . 

[67] M. Isasa, E. Villamor, L. E. Hueso, M. Gradhand, and 

F. Casanova, Physical Review B 91, 024402 (2014), 
arXiv: 1407.4770 . 

[68] E. Saitoh, M. Ueda, H. Miyajima, and G. Tatara, Ap¬ 
plied Physics Letters 88, 182509 (2006). 

[69] K. Ando, S. Takahashi, K. Harii, K. Sasage, J. Ieda, 
S. Maekawa, and E. Saitoh, Physical Review Letters 
101, 036601 (2008). 

[70] T. Kimura, Y. Otani, T. Sato, S. Takahashi, and 
S. Maekawa, Physical Review Letters 98, 156601 (2007), 
arXiv:cond-mat/0609304 . 

[71] K. Uchida, J. Xiao, H. Adachi, J. Ohe, S. Takahashi, 
J. Ieda, T. Ota, Y. Kajiwara, H. Umezawa, H. Kawai, 

G. E. W. Bauer, S. Maekawa, and E. Saitoh, Nature 
materials 9, 894 (2010), arXiv: 1009.5766 . 


18 


[72] Y. Kajiwara, K. Harii, S. Takahashi, J. Ohe, K. Uchida, 
M. Mizuguchi, H. Umezawa, H. Kawai, K. Ando, 
K. Takanashi, S. Maekawa, and E. Saitoh, Nature 464, 
262 (2010). 

[73] A. Aronov and Y. Lyanda-Geller, JETP Letters 50, 431 
(1989). 

[74] V. Edelstein, Solid State Communications 73, 233 
(1990). 

[75] A. Y. Silov, P. A. Blajnov, J. H. Wolter, R. Hey, K. H. 
Ploog, and N. S. Averkiev, Applied Physics Letters 85, 
5929 (2004). 

[76] K. Shen, G. Vignale, and R. Raimondi, Physical Review 
Letters 112, 096601 (2014), arXiv:1311.6516 . 

[77] S. Ganichev, E. Ivchenko, V. Bel’kov, S. Tarasenko, 
M. Sollinger, D. Weiss, W. Wegscheider, and W. Prettl, 
Nature 417, 153 (2002). 

[78] J. C. Rojas Sanchez, L. Vila, G. Desfonds, S. Gam- 
barelli, J. P. Attane, J. M. De Teresa, C. Magen, and 
A. Fert, Nature communications 4, 2944 (2013). 

[79] V. M. Edelstein, Physical Review Letters 75, 2004 
(1995). 

[80] V. M. Edelstein, Physical Review B 72, 1 (2005). 

[81] S. K. Yip, Physical Review B 65, 144508 (2002), 
arXiv:cond-mat/0110140 . 

[82] E. Bauer and M. Sigrist, Non-Centro symmetric Su¬ 
perconductors , Lecture Notes in Physics, Vol. 847 
(Springer, Berlin, Heidelberg, 2012). 

[83] I. V. Krive, A. M. Kadigrobov, R. I. Shekhter, and 
M. Jonson, Physical Review B 71, 214516 (2004), 
arXiv:cond-mat/0409063 . 

[84] A. I. Buzdin, Physical Review Letters 101, 107005 
(2008), arXiv:0808.0299 . 

[85] J.-F. Liu and K. Chan, Physical Review B 82, 125305 

( 2010 ). 

[86] A. G. Mal’shukov, S. Sadjina, and A. Brataas, Physical 
Review B 81, 060502 (2010). 

[87] A. A. Reynoso, G. Usaj, C. A. Balseiro, D. Feinberg, 
and M. Avignon, Physical Review B 86, 214519 (2012), 
arXiv:1212.2786 . 

[88] T. Yokoyama, M. Eto, and Y. V. Nazarov, Journal 
of the Physical Society of Japan 82, 054703 (2013), 
arXiv:1212.5390 . 

[89] T. Yokoyama, M. Eto, and Y. V. Nazarov, (2014), 
arXiv: 1402.0305 . 

[90] T. Yokoyama and Y. V. Nazarov, Europhysics Letters 
(EPL) 108, 47009 (2014), arXiv: 1408.4034 . 

[91] G. Campagnano, P. Lucignano, D. Giuliano, and 
A. Tagliacozzo, (2014), arXiv: 1408.3953 . 

[92] F. S. Bergeret and I. V. Tokatly, Europhysics Letters 
(EPL) 110, 57005 (2015), arXiv: 1409.4563 . 

[93] F. Konschelle, (2014), arXiv: 1408.4533 . 

[94] I. Kulagina and J. Linder, (2014), arXiv: 1406.7016 . 

[95] S. V. Mironov, A. S. Mel’nikov, and A. I. Buzdin, 
(2014), arXiv: 1411.1626 . 

[96] V. Braude and Y. V. Nazarov, Physical Review Letters 
98, 077003 (2007), arXiv:cond-mat/0610037 . 

[97] Y. Tanaka, A. Golubov, S. Kashiwaya, and M. Ueda, 
Physical Review Letters 99, 037005 (2007). 

[98] M. Eschrig, T. Lofwander, T. Champel, J. C. Cuevas, 
J. Kopu, and G. Schon, Journal of Low Temperature 
Physics 147, 457 (2007), arXiv:cond-mat/0610212 . 

[99] R. Grein, M. Eschrig, G. Metalidis, and G. Schon, 

Physical Review Letters 102, 227005 (2009), 

arXiv:0904.0149 . 


[100] J.-F. Liu and K. Chan, Physical Review B 82, 184533 

( 2010 ). 

[101] I. Margaris, V. Paltoglou, and N. Flytzanis, Journal 
of physics. Condensed matter : an Institute of Physics 
journal 22 , 445701 (2010). 

[102] Y. Tanaka, T. Yokoyama, and N. Nagaosa, Physical 
Review Letters 103, 107002 (2009), arXiv:0907.2088 . 

[103] F. Dolcini, M. Houzet, and J. S. Meyer, Physical Re¬ 
view B 92, 035428 (2015), arXiv:1503.01949 . 

[104] V. Geshkenbeim and A. I. Larkin, JETP Letters 43, 395 
(1986). 

[105] S. Yip, Physical Review B 52, 3087 (1995). 

[106] Y. Tanaka and S. Kashiwaya, Physical Review B 56, 
892 (1997). 

[107] M. Sigrist, Progress of Theoretical Physics 99, 899 
(1998). 

[108] S. Kashiwaya and Y. Tanaka, Reports on Progress in 
Physics 63, 1641 (2000). 

[109] Y. Asano, Y. Tanaka, M. Sigrist, and S. Kashiwaya, 
Physical Review B 67, 184505 (2003). 

[110] Y. Asano, Y. Tanaka, M. Sigrist, and S. Kashiwaya, 
Physical Review B 71, 214501 (2005). 

[111] P. Brydon, B. Kastening, D. Morr, and D. Manske, 
Physical Review B 77, 104504 (2008), arXiv:0709.2918 

[112] Y. Rahnavard, D. Manske, and G. Annunziata, Physi¬ 
cal Review B 89, 214501 (2014). 

[113] A. A. Reynoso, G. Usaj, C. Balseiro, D. Feinberg, 
and M. Avignon, Physical Review Letters 101, 107001 
(2008), arXiv:0808.1516 . 

[114] A. Zazunov, R. Egger, T. Martin, and T. Jonck- 
heere, Physical Review Letters 103, 147004 (2009), 
arXiv:0909.3036 . 

[115] A. Brunetti, A. Zazunov, A. Kundu, and R. Egger, 
Physical Review B 88, 144515 (2013), arXiv:1305.3816 

[116] A. I. Buzdin and A. Koshelev, Physical Review B 67, 
220504 (2003). 

[117] E. Goldobin, D. Koelle, R. Kleiner, and A. I. Buzdin, 
Physical Review B 76, 224523 (2007). 

[118] E. Goldobin, D. Koelle, R. Kleiner, and R. G. 
Mints, Physical Review Letters 107, 227001 (2011), 
arXiv:1110.2326 . 

[119] H. Sickinger, A. Lipman, M. Weides, R. G. 
Mints, H. Kohlstedt, D. Koelle, R. Kleiner, and 
E. Goldobin, Physical Review Letters 109, 107002 
(2012), arXiv: 1207.3013 . 

[120] R. Mints, Physical Review B 57, R3221 (1998). 

[121] A. I. Buzdin, Journal de Physique IV (Proceedings) 
131, 213 (2005). 

[122] A. I. Buzdin, Physical Review B 72, 100501 (2005). 

[123] F. Konschelle, The European Physical Journal B 87, 
119 (2014), arXiv:1403.1797 . 

[124] E. Rashba, Physical Review B 68, 241315 (2003), 
arXiv:cond-mat/0311110 . 

[125] I. V. Tokatly, Physical Review Letters 101, 106601 
(2008), arXiv:0802.1350 . 

[126] D. Agterberg and R. Kaur, Physical Review B 75, 
064511 (2007), arXiv:cond-mat/0612216 . 

[127] A. G. Mal’shukov, L. Y. Wang, C. S. Chu, and 
K. A. Chao, Physical Review Letters 95, 146601 (2005), 
arXiv:cond-mat/0506724 . 

[128] T. D. Stanescu and V. Galitski, Physical Review B 75, 
125307 (2007), arXiv:cond-mat/0611165 . 


19 


[129] M. Duckheim, D. L. Maslov, and D. Loss, Physical 
Review B 80, 235327 (2009), arXiv:0909.1892 . 

[130] J. Frohlich and U. Studer, Reviews of Modern Physics 
65, 733 (1993). 

[131] B. Berche and E. Medina, European Journal of Physics 
34, 161 (2013), arXiv:1210.0105 . 

[132] A. Golubov, M. Kupriyanov, and E. Il’ichev, Reviews 
of Modern Physics 76, 411 (2004). 

[133] V. P. Mineev and K. V. Samokhin, Physical Review B 
78, 7 (2008), arXiv:0807.3021 . 

[134] M. Houzet and J. S. Meyer, arXiv (2015), 
arXiv: 1502.04504 . 

[135] G. Eilenberger, Zeitschrift fiir Physik 214, 195 (1968). 

[136] A. I. Larkin and Y. N. Ovchinnikov, Sov. Phys. JETP 
28, 1200 (1969). 

[137] K. D. Usadel, Physical Review Letters 25, 507 (1970). 

[138] N. B. Kopnin, Theory of nonequilibrium superconductiv¬ 
ity (Oxford University Press, 2001). 

[139] C. Gorini, P. Schwab, R. Raimondi, and A. L. 
Shelankov, Physical Review B 82, 195316 (2010), 
arXiv: 1003.5763 . 

[140] D. N. Langenberg and A. I. Larkin, Nonequilibrium su¬ 
perconductivity (North-Holland, 1986). 

[141] A. Mal’shukov and C. Chu, Physical Review B 78, 
104503 (2008), arXiv:0801.4419 . 

[142] M. Y. Kupriyanov and V. F. Lukichev, Sov. Phys. JETP 
67, 1163 (1988). 

[143] L. P. Gor’kov and E. Rashba, Physical Review Letters 
87, 37004 (2001). 

[144] I. Sochnikov, L. Maier, C. A. Watson, J. R. Kirt- 
ley, C. Gould, G. Tkachov, E. M. Hankie wicz, 


C. Briine, H. Buhmann, L. W. Molenkamp, and K. A. 
Moler, Physical Review Letters 114, 066801 (2014), 
arXiv:1410.1111 . 

[145] A. Bauer, J. Bentner, M. Aprili, M. L. Della Rocca, 
M. Reinwald, W. Wegscheider, and C. Strunk, Physical 
Review Letters 92, 4 (2003), arXiv:cond-mat/0312165 . 

[146] A. Buzdin, Journal of Physics: Conference Series 150, 
052029 (2009). 

[147] L. Serrier-Garcia, J. C. Cuevas, T. Cren, C. Brun, 
V. Cherkez, F. Debontridder, D. Fokin, F. S. Bergeret, 
and D. Roditchev, Physical Review Letters 110, 157003 
(2013), arXiv:1401.8102 . 

[148] F. Konschelle and A. I. Buzdin, Physical Review Letters 
102, 017001 (2009), arXiv:0810.4286 . 

[149] Lower indices (i, j, k...) describe space variables, while 
upper indices (a, 5, ...) will describe spin variables. 

[150] We use the units where the Planck’s constant h— 1 and 
the Boltzmann’s constant ks — 1. 

[151] Note that, in a one dimensional system described by 
the Hamiltonian Eq. (3.2) the gauge-potential is al¬ 
ways a pure-gauge and therefore ID problem should 
not exhibit any magneto-electric effect. Nevertheless, 
the above argument does not apply for topological ID 
channel, where a (^o-effect has been discussed [102, 103], 
since the dispersion relation is not quadratic in that 
case [95]. Also, the problem of spin-active interface is 
non-trivial when spin-orbit effect is treated as a gauge- 
potential, and might eventually lead to (^o-effect as well 
[94, 99-101]. 


